Exotic Kondo effect from magnetic trimers 
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Motivated by the recent experiments of Jamneala et al. [T. Jamnaela et al., Phys. Rev. Lett. 87, 
256804 (2001)], by combining ab-initio and renormalization group methods, we study the strongly 
correlated state of a Cr trimer deposited on gold. Internal orbital fluctuations of the trimer lead 
to huge increase of Tk compared to the single ion Kondo temperature explaining the experimental 
observation of a zero-bias anomaly for the trimers. The strongly correlated state seems to belong 
to a new, yet hardly explored class of non-Fermi liquid fixed points. 
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PACS numbers: 75.20.Hr, 71.27.+a, 72.15.Qm 

In recent years, atomic scale resolution Scanning Tun- 
neling Microscopy (STM) proved to be a spectacular tool 
to probe the local density of states around Kondo impu- 
rities adsorbed on a metallic surface 0, . Experiments 
were performed with single Ce atoms on Ag and Co 
atoms on Au and Cu surfaces 0,0]. When the STM 
tip is placed directly on the top of the magnetic adatom, 
a sharp resonance appears in the differential conductance 
at low bias and disappears when the STM tip is moved 
away from the impurity or when the substrate temper- 
ature is raised above the Kondo temperature Tk- This 
zero bias anomaly appears as the main signature of the 
Kondo effect and results from the screening of the adatom 
spin by the surrounding (bulk and surface) conduction 
electrons. The precise line shape of the resonance can 
well be understood in terms of a Fano resonance, an in- 
terference phenomenon occurring because of two possible 
tunneling channels: a direct channel between the tip and 
the impurityand a second channel between the tip and 
the surface 0,11 00. 

The manipulation of single atoms on top of a surface 
with an STM tip has also been proven useful to build clus- 
ters of atoms with well-controlled interatomic distances. 
For example, Manoharan et al. manufactured an el- 
liptical quantum corral of Co adatoms, and found that, 
when an extra Co adatom is placed at one focus of the 
elliptical corral, a "mirage" of the Kondo resonance can 
also be observed at the other focus. Another intriguing 
result was found recently for Cr trimers on a gold surface 
by Jamneala et al. 0. Whereas isolated Cr monomers 
or dimers display featureless signals in STM spectra at 
T = 7K, Cr trimers exhibit two distinct electronic states 
depending on the atomic positions: a sharp Kondo reso- 
nance of width Tk ~ 50K was found for an equilateral 
triangle, while the STM signal of an isosceles triangle 
did not show any particular feature. Furthermore, the 
Cr trimers were reversibly switched from one state to 
another. As schematically depicted in Fig. ^ Cr atoms 
forming an equilateral triangle are expected to occupy 



nearest neighbor sites on the gold surface (Au(lll)), al- 
lowing therefore geometric frustration |Sj. Such a com- 
pact magnetic nanocluster is of very much theoretical 
interest due to the interplay between Kondo physics and 
magnetic frustration that may generate internal orbital 
fluctuations. This system can indeed be regarded as the 
smallest and simplest frustrated Kondo lattice and, as 
we will see, embodies a very rich behavior. We show in 
this Letter that internal orbital fluctuations in the Cr 
equilateral trimer lead to a huge increase of the Kondo 
temperature in agreement with the experiment and we 
argue that the low energy physics of the system is gov- 
erned by a new non-Fermi liquid fixed point 0. 

In this Letter we shall present a careful study of the 
strongly correlated state of the Cr trimer depicted in 
Fig. ^ I n order to construct an effective Hamiltonian 
for this system, we first performed ab-initio calculations 
\vl\ to study the electronic structure of the Cr ions form- 
ing the cluster, and verified that the Cr ions are within 
a very good approximation in d 5 spin Scr = S = 5/2 
states, and display relatively small valence fluctuations. 
Under these conditions, the Hamiltonian describing the 
nanocluster can well be approximated as 
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H = H spin + — S l i\)\ (J a (J(jl vb i , 



(1) 



where H sp i n describes the interaction between the Cr 
spins Si at sites i = 1,..,3, and G denotes the Kondo 
coupling between each Cr spin and the conduction elec- 
trons in the substrate. Note that Eq. Q incorporates 
only the exchange generated by the most strongly hy- 
bridizing d-state. This approximation is justified by 
the observation that the Kondo effect is exponentially 
sensitive to the strength of hybridization, and there- 
fore hybridization with other <i-states is not expected 
to change our results. Correspondingly, it is sufficient 
to consider only the hybridization with a single conduc- 
tion electron state, ipi, leading to Eq. JJJ. Note, how- 
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a rather simple form in this notation, 
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FIG. 1: Fig. a: Top view of the equilateral Cr cluster on 
the Au(lll) surface. Au and Cr atoms are indicated by grey 
and black circles, respectively. Fig. b: Eigenstate structure 
of //spin of the Cr cluster (1) in the absence and (2) in the 
presence of spin-orbit coupling. 



ever, that these states may overlap with each-other and 
in general {ipi^ip]} + 1 As we demonstrate later by 
ab-initio calculations, the effective spin-spin interaction 
ffspin in Eq. Q is rather well approximated by a dipolar 
term ff spin w ffdipoi = JJ2(i,j) where the coupling 

J w 1600K turns out to be antiferromagnetic. 

Since the Kondo temperature Tk ~ 50K generated 
by G is expected to be much smaller than the exchange 
coupling J, and we only wish to describe the physical be- 
havior of the cluster around the low energy scale ~ Tk, 
we shall first diagonalize ff sp i n and construct an effec- 
tive Hamiltonian to describe quantum fluctuations of the 
cluster spins. The The low energy section of the spectrum 
of -ffdipoi is sketched in Fig. ^p: There are two different 
ways to construct states with total spin S — 1/2 that 
minimize -ffdipoi- As a consequence, the ground state is 
fourfold degenerate, the extra degeneracy being associ- 
ated with a two-dimensional representation of the C$ v 
symmetry of the cluster. These four spin states can thus 
be labeled as |<r, //), where a =t, j denotes the spin and 
H = ± is a chiral index: Ca\a,fJ,) = e v27r / 3 |cr, fi). Note 
that the four states remain degenerate as long as ff sp i n 
is SU(2)-invariant, and only spin-orbit coupling related 
effects discussed later can split their degeneracy. 

Next , we shall obtain the effective Hamiltonian within 
the subspace {\a,fi}}. To do this we first construct the 
states | cr, n) explicitly using Clebsch-Gordan coefficients 
and then evaluate the matrix elements of the Kondo ex- 
change, H K = yX)<,tr,o-' Siipi^a^/tpia' within this sub- 
space in order to obtain an effective Hamiltonian ff e ff. 
The effect of virtual transitions into the high energy sub- 
space on the effective low energy couplings is neglected 
assuming G/J-^i 1 To express the effective Hamiltonian 
ffeff in a simple form we define the following Fcrmionic 
fields, tpvp = ^^j e ' w 2f ' 3 ^ji" where now fi takes 
three possible values, /i = 0, ±. After tedious and lengthy 
algebraic manipulations, the effective Hamiltonian takes 
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- (25+ l)(ST+ ip^T-ip + h.c. 



(2) 



where S denotes a spin 1/2 operator acting on the spin 
indices of \<j,/j,), and the orbital pseudospin operators 
T ± are standard operators raising/lowering the chiral 
spin \i. The operators r in Eq. @ change the angu- 
lar quasi momentum of the conduction electrons by one 
unit, r^„/ = tfffL-tii where 6^ denotes the Kronecker 
delta function modulo 3. Note that ff e // depends on the 
magnitude of the initial local spins enhancing the second 
term by a factor 25+1. 

To determine the low-energy dynamics of the system 
and the Kondo scale Tk associated with the binding en- 
ergy of its strongly correlated state, we carried out a 
perturbative renormalization group (RG) analysis of the 
system. Within the RG we integrate out conduction elec- 
trons with large energy and take their effect into account 
by changing the coupling constants of the original model. 
Thereby, the original bandwidth Dq is gradually reduced 
to smaller and smaller values D. Two stages must be dis- 
tinguished in course of the RG procedure: (a) For D 3> J 
the neighboring Cr spins behave as independent spins, 
(b) For D <C J on the other hand the Cr spins are tied 
together by their exchange interaction, and the effective 
Hamiltonian can be used. In the first regime only G is 
renormalized according to the usual scaling equation^] : 



dG - n r 2 1 n 2 r 3 
~df- goG 2 0oG 



(3) 



where I = \n(Do/D) is the scaling parameter and go de- 
notes the local density of states, and is related to the 
propagator of the field tp a j through (ip a j ;(t)ip} r j (0)) f» 
-igo/t. 

At energies (time scales) below J, we can thus use the 
effective Hamiltonian but with a conduction elec- 
tron cut-off reduced to D — D ~ 3 J/2 and the coupling 
G replaced by an effective coupling G — > G = G(l — 
\n(Do/Do)). In this regime, however, the RG analysis 
becomes more complicated, and additional terms are gen- 
erated in the Hamiltonian. To allow for the generation of 
these terms we first introduce a very general Hamiltonian 
of the form: 



(4) 



a,/3=l,..,4p,g=l,..6 



where a, = {a = ±,fi = +} and p, q = {a = 
±, H = ±,0} denote composite indices referring to the 
ground state multiplet and the conduction electrons, re- 
spectively. To obtain Tk we solved numerically the RG 
equations derived in Ref. [l3| for this general model. 
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FIG. 2: Kondo temperature as a function of the overlap 
parameter a for a cluster formed of atoms with spin S = 
5/2 and spin S = 1/2. We used D = 5000 K « 3 J and a 
dimensionless coupling G = 0.09. The inset shows in a log- 
log plot the Kondo temperature of the trimer as a function of 
the single Cr's Kondo temperature for two overlap paraemetrs 
a = 0.33 and a — 0.6. 



The initial values of the couplings V^f can easily be 
determined from Eq. J3J) at D = Dq. However, spe- 
cial care is needed to define the dimensionless couplings 
entering the scaling equations: since {tpi,ipj} ^ <5,.j, 
the off-diagonal-correlation function decays in general as 
("i/viWV^W) ~ ~i a Qo/t> with an overlap parameter 
a < 1. As a result, the density of states in the electronic 
channels depends on the chiral index: g^ = Qo(l + 2a) 
for fi — and g^ — go(l — a) for /i = ±. These density 
of states enter the dimensionless couplings of Ref. [13| as 
v pq = y/QpQq Vp<f ■ The parameter a can be estimated 
using a free electron model [lfj , but a simple tight bind- 
ing model can also be used to estimate it |l9j], and is 
typically in the range 0.2 < a < 0.6. 

The advantage of the RG method discussed above is 
that it sums up systematically leading and subleading 
logarithmic corrections, and provides an unbiased tool to 
obtain Tk, even without knowing the structure of the 
fixed point Hamiltonian. This is in sharp contrast with a 
variational study as the one developed by Kudasov and 
Uzdin for this problem |l4|. which only gives an upper 
bound on Tk, gets even the exponent of Kondo temper- 
ature incorrectly (see e.g. in Ref. an d the result 
of which depends essentially on the specific variational 
ansatz made. Furthermore, the perturbative renormal- 
ization group enables one to gain some insight into the 
structure and symmetry of the fixed point Hamiltonian 
describing the physics of Eq. (JTJ below the Kondo scale. 

The Kondo temperature can be defined as the en- 
ergy scale at which the norm of the effective couplings, 
N = (Eo 8 p q \ v pq I 2 ) 1 / 2 reaches the strong coupling 
limit (N{D = T K ) = 0.7). For a = 1 only the 
field Vv=o couples to the cluster, and the scaling equa- 
tions reduce to those of an isolated Cr spin. Therefore, 



in this limit Tk is the same as that of an isolated Cr 
ion, T^\ In the inset of Fig. H we show the Kondo 
temperature of the trimer as a function of the Kondo 
temperature of a single Cr ion for two typical values of 
the overlap parameter a and for an intermediate cut-off 
Dq = 5000AT ~ 3J. The trimer's Kondo temperature 
is orders of magnitude larger than that of the single Cr 
ion, and for 0.001-tT < Tjp < OAK (which is in agree- 
ment with the small bulk Kondo temperature of Cr in Au 
[III), yields a trimer Kondo scale 10K < T K < WOK. 

We also find that T K « C(a) \J Dq T^ , with C(a) a 
constant of the order of unity. This increase is clearly 
due to the presence of orbital fluctuations, and gives a 
natural explanation to the fact that the Kondo effect has 
not been observed experimentally for a single Cr ion on 
Au.|j| Note, however, that this increase is peculiar to 
the case of large cluster spins S, and no such a dramatic 
increase appears for spin 5 = 1/2 trimers, in contrast 
to the results of Ref. [1J|. This important difference of 
behavior between S — 1/2 and S — 5/2 is highlighted 
in Fig. |2 where on a log scale we displayed Tk as a 
function of the overlap parameter for the specific choice, 
G = 0.09, corresponding to a trimer Kondo scale ~ 40 K 
and re 0.01 K. 

Distorting the equilateral trimer configuration, we 
clearly lift the 4-fold ground state degeneracy of -ffspin 
and therefore suppress this Kondo cluster effect and no 
enhancement of the Kondo temperature is expected. Ge- 
ometric distortion in our formalism appears in fact as a 
strong orbital magnetic field. This gives a natural expla- 
nation why the isosceles Cr trimer displayed no Kondo 
resonance at the experimental temperature T = 7K . 

The RG flows also allow us to gain information on 
the symmetry of the low energy fixed point, governing 
the T — > dynamics of the cluster, since dynamically 
generated symmetries typically show up already within 
the perturbative RG. Surprisingly enough, analyzing the 
structure and algebraic properties of the fixed point cou- 
plings , we find that it is neither the familiar SU(4) 
Fermi liquid fixed point [lij nor the two-channel Kondo 
fixed point ^t|- We find instead that the fixed point 
Hamiltonian takes the following from: 

Hfp = a Sa + (3 Sa L 2 Z +-/T Z L Z 

+ S Sa (T-L+ +T+L-) , (5) 

where for simplicity we suppressed the fermion fields. 
The operators L ± and L z above denote standard L = 1 
angular momentum matrices acting on the fermionic or- 
bital spin. While this Hamiltonian clearly has an en- 
larged U(l) symmetry [replacing the original C3 symme- 
try], it remains nevertheless anisotropic in orbital space. 
This structure allows us to identify this fixed point with 
the non- Fermi liquid fixed point identified first in Ref. 
using the numerical renormalization group method for 
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FIG. 3: Energy of the Cr trimer as a function of the direc- 
tion of magnetization computed using a relativistic ab-initio 
method (full line) and the effective model, Eq. (dotted 
line) for the particular configurations sketched in the inset. 
In this example only the angle & with respect to the z axis is 
varied (upper inset). The lower inset shows the projections 
of the magnetization onto the surface (xy plane) . 

the simpler problem of a spin 1/2 model trimer. A sup- 
plementary strong coupling analysis can be used to sup- 
port that this fixed point is indeed at finite couplings and 
must therefore be of non-Fermi liquid character. We also 
find that this fixed point is unstable to general hermi- 
tian perturbations (another indication of unstable fixed 
points), under which it flows to the familiar stable SU(4) 
Fermi liquid fixed point [15j. On the other hand, we could 
not find additional terms respecting SU (2) x C% v symme- 
try that would render the strange fixed point unstable. 
The analysis of the physical properties of this non-Fermi 
liquid fixed point is beyond the scope of the present pa- 
per, but the thermodynamical, transport, and tunneling 
(STM) properties are expected to show anomalous scal- 
ing properties. 0,0] 

So far we fully neglected the effect of spin-orbit cou- 
pling. However, spin-orbit coupling removes the four-fold 
degeneracy of the ground state of H sp - mi and splits it up 
into two Kramers doublets. To build a consistent picture 
(and also to have a non-Fermi liquid regime) this split- 
ting A must be smaller than Tk- In order to estimate 
the splitting we performed relativistic ab-initio calcula- 
tions ^3 by fixing the orientation of Cr spins, Sj — > SQj, 
and determined the energy of the cluster as a function of 
the spin orientations, i?[{f2j}]. Details of this calcula- 
tions will be published elsewhere. 0] The energy of the 
cluster can excellently be approximated by the following 
expression, 

= S 2 ]T Jf^ + H Q . (6) 

a,f3—x,y,z 

The dominant term of the interaction turns out to be the 
simple SU(2) invariant exchange interaction studied in 



the first part of the paper. The term Hq above contains 
SU(2) invariant quadrupolar couplings and three-spin in- 
teractions. While these turn out to be of the same order 
of magnitude as the anisotropy of the exchange interac- 
tion, Hq does not lead to a splitting of the ground state 
degeneracy in contrast to anisotropy terms. In Fig. E| the 
ab-initio values of the energy of the cluster are compared 
with energies of the effective Hamiltonian in Eq. @ for 
a particular configuration of the Cr spins. Having de- 
termined the exchange couplings above, we performed 
first order perturbation theory within the degenerate sub- 
space of the trimer to obtain A ss 20K, which is indeed 
less than the experimentally observed Kondo tempera- 
ture, Tk ~ 50K. This splitting will ultimately destroy 
the non-Fermi liquid properties predicted above, how- 
ever, the energy scale A* at which this happens is ex- 
pected to be smaller than A similar to the two-channel 
Kondo model, where A* ~ A 2 /T K . 

For Cr on gold the ratio TV /A* is probably too small 
to observe the non-Fermi liquid physics. However, it 
should be possible to observe it in other systems: Cr on 
Ag, e.g., is a promising candidate since the lattice con- 
stant of Ag is about the same as that of Au (important 
to have in antiferromagnetic coupling between magnetic 
ions), while the spin-orbit coupling is much weaker. Our 
calculations show that in this case A ~ IK, and a much 
wider non-Fermi liquid range may be accessible. In this 
case the atoms would be replaced by quantum dots. The 
great advantage of such a device would be that (1) it 
would be highly tunable and (2) it would allow for trans- 
port measurements, though it is not easy to guarantee 
the perfect symmetry of the device. 
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